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1 Introduction 



With the fruitful running of both BaBar and Belle in the past decade, as well as the upcoming 
LHC-b and the proposed Super-B experiments, rare B-meson decays provide a golden oppor- 
tunity to test the Standard Model (SM) picture of flavor physics and CP violation. Although 
most experimental measurements are in perfect agreement with the SM predictions, some in- 
comprehensible discrepancies, such as the measured large branching ratio of B° — > 7r 7r° [I] (the 
so-called "irir puzzle") and the direct CP asymmetries Acp{B~ — > ir K~) ^ Acp(B° — > n + K~) 
at 5cr significance [2] (the so-called "ttK puzzle" OH]), still exist. 

The observed u inr puzzle" is reflected by the following two ratios of CP-averaged branching 
fractions [HJ 0] 



B(B~ ->7r-7r°) -ir B o r B(B° ttV 

B(B° -)■ 7T+TT-)\ T B+ ' °° LB(B° ->■ 7T+7T- 



(1) 



With the up-to-date results averaged by the Heavy Flavor Averaging Group (HFAG) [T] from 
BaBar Belle [SI E], CLEO [TO] and CDF [TJ], and taking r B +/r B o = 1.071 ±0.009 [12], 

we get 

R%L (Exp.) = 2.02 ±0.17, ^(Exp.) = 0.60 ±0.08. (2) 

Theoretically, it is generally expected that B(B° -»■ 7r + 7r _ ) > Z3(£?~ ->■ 7r~7r°) and B(B° -)• 
7r + 7r _ ) ^> B(S° — > 7r°7T ) within the SM. Adopting the central values calculated within the 
QCD factorization (QCDF) pH El US] , the perturbative QCD (pQCD) [TBI El US] and the 
soft-collinear theory (SCET) [T9J [2D1 EEEJ [22] , we get respectively 

#™ (QCDF) = 1.25, 1.83 (S4) , i^(QCDF) = 0.07, 0.27 (S4) [23] , 

fl™ (pQCD) = 0.93, 1.15 (+NLO) , i^(pQCD) = 0.03, 0.09 (+NLO) [21] , (3) 
iCTL(SCET) = 1.80, i^(SCET) = 0.31 



Comparing with the experimental data Eq. ([2]), we can see that (i) within the Scenario S4 of 
QCDF and the SCET formalism, even though the ratio R™_ is well predicted, their predictions 
for Rqq are still much lower than the data; (ii) the pQCD predictions are, on the other hand, 
even worse no matter the NLO corrections are included or not. It is also noted that, although a 
large electro-weak penguin (EWP) contribution could resolve the anomaly observed in B — > tttt 



and ttK decays [31 HJ HHJ [27] , there are no known mechanisms to enhance such a large EWP 
amplitude. Furthermore, a large color-suppressed tree amplitude, which is also helpful to 
moderate these puzzles [2H EH [291 EDI EI] , is hard to obtain from any short- distance dynamics. 
Thus, it is quite difficult to understand these anomalies within the SM. Many efforts have been 
made to bridge these large discrepancies both within the SM [321 EH EU ESI ESJ EJJ EH EH ED] 
and in various new physics (NP) scenarios [%T] H21 H31 HH H5] . 

As is well-known, an additional U(l)' gauge symmetry and associated Z' gauge boson could 
arise in some well-motivated extensions of the SM. Searching for such an extra Z' boson is an 
important mission in the experimental programs of Tevatron [16] and LHC (UJ HH] . Performing 
the constraints on the new Z' couplings through low-energy physics is, on the other hand, 
very crucial and complementary for these direct searches. Theoretically, one of the simple 
extensions is the family non-universal Z' model, which could be naturally derived in certain 
string constructions (121 [501 EE] , E% models [521 [531 EH E21 [56] and so on. It is interesting to note 
that the non-universal Z' couplings could lead to large flavour- changing neutral current (FCNC) 
processes and possible new source of CP- violating effects [571 EH] - Such a specific model could 
reconcile the anomalies observed in B s — B s mixing and B — > <pK s decay [591 EDI EH E2] , the 
li ttK puzzle" [631 16^1 165] , and the anomalous polarization in B — > <pK* decays [66] . It has also 
been examined in some other interesting processes [63, EB1 ESI EDI EH [721 [731 EU ES] • 

Based on the above observations, in this paper we pursue a possible solution to the ob- 
served "7T7r" puzzle within such a family non-universal Z' model. We shall adopt the QCDF ap- 
proach [131 [HI [T51 E3] to evaluate the relevant hadronic matrix elements of B — > tttt decays, with 
an alternative scheme to parameterize the end-point divergence appearing in hard-spectator 
and annihilation corrections [6H [76] . In addition, since only the flavour- nondiagonal couplings 
B P b (p = d, s) are different between the quark-level transitions b — > dqq and b — > sqq (with 
q = u or d) within the model, it is interesting to investigate if the allowed parameter spaces 
in B — > tttt decays could survive after taking into account the constraints from B — > irK, ttK* 
and pK decays [61]. 

Our paper is organized as follows. In Section 2, a survey of B — > tttt decays in the SM within 
the QCDF formalism is given; our numerical result, with two different schemes for the end-point 
divergence, is also presented. In Section 3, after briefly reviewing the non-universal Z' model, 
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we present our numerical results and discussions in detail. Section 4 contains our conclusions. 
Appendix A recapitulates the SM decay amplitudes of the B — > tttt decays within the QCDF 
formalism [131 H2 EH E3] . Appendix B contains the relevant formulas for hard-spectator and 
annihilation amplitudes with the infrared finite gluon propagator [BU [TB] . All the theoretical 
input parameters are summarized in Appendix C. 



2 Revisiting B — > tttt decays within the QCDF framework 

In the SM, the effective weak Hamiltonian responsible for b — > d transition is given as [771 EH] 



n/SM _ 



r r 10 
v ub v: d (CiO? + c 2 o u 2 ) + v cb v; d {doi + c 2 o 2 c ) - v tb v; d ( £ ao, 



i=3 



CjryOjry + C $gO gg) 



-h.c, (4) 



where V q bV* d (q = u, c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements [7H1 ED] , Q the Wilson coefficients, and Oj the relevant four-quark operators whose 
explicit forms could be found, for example, in Refs. [131 El US Ell ITS] . 

In recent years, the QCDF approach has been employed extensively to study hadronic 
B-meson decays. For example, the tree-dominated B tttt decays have been studied compre- 
hensively within the SM in Refs. [131 EH 123] . It is also noted that the framework contains 
estimates of the hard-spectator and annihilation corrections. Even though they are power- 
suppressed, their strengthes and associated strong-interaction phases are numerically impor- 
tant to reproduce the experimental data [231 ESI 123 EOl EI] • Unfortunately, there are end-point 
divergences associated with the twist-3 hard-spectator and the annihilation amplitudes. Thus, 
how to regulate the divergence becomes indispensable within the QCDF formalism. 

As the usually adopted scheme (Scheme I), the divergent integrals are phenomenologically 
parameterized by X A (annihilation) and Xh (hard spectator scattering) [131 CHI US [23] , with 



^X A = (l + p A e^ A ) In / dy^ --(X A ) 2 , (5) 

o y A h Jo y 2 

and similarly for Xjj . Here A^ = 0.5 GeV, Ph,a < 1 and 4>h,a is an arbitrary strong-interaction 
phase, which might be caused by soft rescattering. The different choices of p A and <fr A correspond 
to the different scenarios discussed in Ref. [23]. In our numerical evaluations, we take Ph,a — 1 



Table 1: The CP-averaged branching ratios (in units of 10 6 ) and the CP asymmetries (in unit of 
10 -2 ) of B — > 7T7T decays in the SM with the two regulation schemes for the end-point divergence. 



Observable 


Exp. data 


Scheme I 


m g = 0.3 GeV 


Scheme II 
m g = 0.7 GeV m g 


= 0.5 ± 0.05 GeV 


B(B~ -)> tt-tt ) 


O-O»_0.40 


4.64 ±0.74 


5.42 ±0.65 


4.01 ±0.64 


4.52 ±0.70 


B(B° 71 + 71-) 


5.16 ±0.22 


6.88 ± 1.13 


9.32 ± 1.07 


7.44 ± 1.00 


8.14 ± 1.21 


B(B° 7T 7T ) 


1.55 ±0.19 


0.77 ±0.12 


2.03 ±0.24 


0.57 ±0.07 


0.95 ±0.18 


A C p(B~ -»■ 7T-7T ) 


6±5 


-0.1 ±0.1 


-0.1 ±0.1 


-0.1 ±0.1 


-0.1 ±0.1 


A CP (5° -»■ 7T + 7T-) 


38 ±6 


15 ±2 


42 ±3 


18 ±2 


26 ±3 


A CP (5° -»• 7r°7r°) 


4 o+25 


-7±6 


-60 ±4 


-35 ±4 


-45 ±5 




-65^ 


-58 ±7 


-69 ±6 


-52 ±6 


-64 ±7 


A£j?(£ -> 7r°7r°) 




70 ± 10 


17 ±4 


45 ±4 


33 ±6 



and 4>h,a = —55°, which are suggested by the most favorable scenario S4 in Ref. [23]. For 
estimating theoretical uncertainties, we shall assign an error of ±0.1 to Ph,a and ±5° to 4>h,a- 
As the above parametrization is rather arbitrary, it is still very worthy to find some al- 
ternative schemes to regulate these endpoint divergences, as precisely as possible, to estimate 
the strengths and the associated strong phases in these power-suppressed contributions. As an 
alternative scheme (Scheme II), we shall use the infrared finite gluon propagator to regulate 
the end-point divergence. It is interesting to note that recent theoretical and phenomenological 
studies are now accumulating supports for a softer infrared behavior of the gluon propaga- 
tor [8TJ[8JI1[S3[S3I8]3I8^ which is also indicated by recent lattice simula- 
tions [921 E3] • Furthermore, the infrared finite dynamical gluon propagator, which is shown to 
be not divergent as fast as 1/q 2 , has been applied successfully to two-body hadronic B-meson 
decays [611 ESI IHH ESI ESI E0 EH]- To be specific, in our calculation we shall adopt the gluon 
propagator derived by Cornwall (in Minkowski space) [221 HOP] 

D{q2) = g 2 -ilW) + , e ' (6) 
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where the dynamical gluon mass square M 2 (q 2 ) is obtained as [221 [100J 



M 2 (q 2 )=m 2 g 



In 



q 2 +im 2 \ -. _ 12 



T 2 



ln<^ 



A 7 

QCD 
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(7) 



with q being the gluon momentum, m g the effective gluon mass and Aqcd = 225 MeV. The 
corresponding strong coupling constant reads 

4-7T 



Pohx 



3 j 

~P 

QCD 



where /?o = H — \fif is the first coefficient of the QCD beta function, and rif the number of 
active quark flavors. 

With such a scheme, the hard-spectator scattering contributions are real, while the anni- 
hilation contributions are complex with a large imaginary part (76J. Both the hard-spectator 
scattering and the annihilation contributions are very sensitive to the value of the effective 
gluon mass pH [98], with a typical value m g = 500 ± 200 MeV [991 EDO]- In Ref - D3], we pre- 
sented our suggested value, m g = 500 ± 50 MeV, which is a reasonable choice so that most of 
the observables for B — > ttK, irK* and pK decays are in good agreement with the experimental 
data. Furthermore, comparing to the available data on B° K+R-^K, B s — > 7T7T and so 
on, Natale and Zanetti also found that the gluon mass scale is close to 500 MeV [98J. So, in 
this paper, we shall take 500 MeV as the central value of the gluon mass scale and ±50 MeV 
as its uncertainty. As a comparison, the numerical results with m g = 700 MeV and 300 MeV 
are also presented in Table [TJ 

With the above two schemes, our SM predictions for the branching ratios, the direct and 
the mixing-induced CP asymmetries of B — > nn decays are listed in Table [TJ We find that, 
although B(B° — > 7r°7r°) in Scheme II could be enhanced relative to the one based on Scheme I, 
B(B° — > 7r + 7r~) is also increased at the same time. Thus, as expected within the SM, we get the 
ratios R 7 ^ = 1.04 ±0.22 and Rqq = 0.23 ±0.06, which are 3.5<r and 3.7<r lower than the current 
experimental data Eq. ([2j), respectively. In addition, our prediction Acp(B° — > 7r°7r°) ~ —0.45, 
although being roughly consistent with the central value in Scheme I ~ —0.07 and that in 
SCET ~ —0.58 [25], is still quite different from the current data 0.431^24 [I]- I n the following, 
we shall investigate if these mismatches could be reconciled within a family non-universal Z' 
model. 
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3 Solution within the family non-universal Z' model 



3.1 Formalism of the family non-universal Z' model 

A family non-universal Z 1 model can lead to FCNC processes even at tree level due to the 
non-diagonal chiral coupling matrix. The formalism of the model has been detailed in Ref. [5T[ 
158] . The relevant studies in the context of B physics have also been extensively performed in 

Refs. [591 EQl EU E21 E31 EH ESI EHl EZ1 ESI [691 [ZQl [TU H21 EH ES]. 

With the assumption of flavor-diagonal right-handed couplings and neglecting the kinetic 
mixing term as usually adopted in the literature, the Z' part of the effective Hamiltonian for b — > 
pqq (p = d,s and q denotes the active quarks) transitions can be written as [57J EH EH EHJ E3] 

n% = ^f(^) 2 B L pb {pb)v-A £ [B L qq {qq)v-A + B*(qq) v+A ] + h.c. , (9) 

where g\ = e/(sin 6w cos 9w), g' is the gauge coupling constant of extra U'{1) group, Mz' the 
mass of the new gauge boson, and refer to the effective Z' couplings to the quarks i and j 
at the electro-weak scale. It is noted that the forms of four-quark operators in Eq. already 
exist in the SM. As a result, we can rewrite Eq. ([9]) as 

n% = -^V tb V t * p J2 + AC 5 O q 5 + AC 7 O q 7 + AC 9 O q 9 ) + h.c. , (10) 

where 0\ [i = 3, 5, 7, 9) are the effective four-quark operators in the SM, and ACi denote the 
modifications to the corresponding SM Wilson coefficients induced by the new gauge boson Z', 
which are expressed as 



ACa 7 



{^) 2 B L pb {BL R -B L dd R ), (11) 



3V tb V* p " gi M z , 

in terms of the model parameters at the initial scale /i^ ~ mw (the W ± boson mass). 

Generally, the diagonal elements of the effective coupling matrices B^' R are real as a con- 
sequence of the hermiticity of the effective weak Hamiltonian. However, the off-diagonal ones 
Bp b can contain a new weak phase (ftp . Then, for convenience we can represent ACi as 



AC 3 , 5 = 2 ^ Ci L ' LR e*# , AC 9J = 4 ^ £™ , (12) 

v tb V tp V tb V tp 



where the newly introduced Z' parameters (p L ' LR , £,p L ' LR and 0^ are defined, respectively, as 



/■ll,lr _ 1 / g'M z n 2 1 B pb . L R ^ L)R . 
Cp " 3 [ g 1 M z/ ) W tb V t ; l[ uu + dd J 



cLL,LR l( 9'Mz \2, B^ b j r-n>L,R f> L ,R\ 

?p 3 { g 1 M z , ) l V tb V t * p l{ uu dd h 

4i = arg[5 p L 6 ]. (13) 

It is noted that the other SM Wilson coefficients may also receive contributions from the Z' 
boson through renormalization group (RG) evolution. With our assumption of no significant RG 
running effect from the scales M' z to Myy, the RG evolution of the modified Wilson coefficients 
is exactly the same as that in the SM [77J EE]- Numerical results of the Wilson coefficients at 



the lower scales m;, and y/A^mb are listed in Appendix CI. 

In our analyses, we also define the following three ratios (taking as a benchmark) 

cLR r>R _ Ri? fLL R L , 9 r>L fLR r>R , o r>R 

r> — D uu D dd d _ D uu ~r z ' D dd p _ D uu "r / ' D dd /-ia\ 

~ f LL ~ B L — B L ' f LL ~ B L — B L ' cll ~ rl _ dl • y iq: ) 

D uu D dd D uu D dd D uu D dd 

It is interesting to note that these ratios are independent of the quark flavor p (p = d, s), i.e., 
they are the same for both b — > sqq and b — > dqq transitions. In fact, there are only two different 
Z' parameters \B p i\ and 0^ between these two types of transitions. Thus, the Z' parameter 
spaces constrained by B — > -kit decays should also suffer constraints from B — > irK, ttK* and 
pK decays, which have been investigated in detail in our previous paper [53] . 

3.2 Numerical results and discussions 

With the theoretical formulae and the input parameters summarized in Appendices A, B and C, 
we now present our numerical analyses and discussions. We take the current experimental data 
on B(B — > titt) and Aqp{B — > tttt) as constraints and A™p(B — > 7r + 7r~, 7r°7r°) as our theoretical 
predictions. Our fitting is performed with the experimental data varying randomly within 2a 
error-bars, while the theoretical uncertainties are obtained by varying the input parameters 
within the regions specified in Appendix C. For simplicity, in the analysis and discussion of this 
subsection, we mainly adopt the Scheme II with m g = 0.5 ± 0.05 GeV to regulate the end-point 
divergence. While, as a comparison, the fitting results with scheme I are also presented. 
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Table 2: Numerical results for the Z' coupling parameters ^ L ' LR , (j R and <f>^ in the three different 
specific cases with the regulation scheme II (scheme I) for the end-point divergence. The dashes in 
Case I means that the corresponding parameters are neglected. 

Parameters Case I Case II Case III 

-5.8 ±6.3 (-5.9 ±4.9) 7.6 ± 1.2 
-13.2 ±7.5 (-13.5 ±8.6) -2.4 ±1.7 
48.5 ± 5.2 (37.0 ± 5.4) 49.5 ± 6.9 (41.8 ± 9.0) 45.4 ± 5.0 
-51 ±5 (-71 ±6) -52 ±6 (-68 ±7) -48 ± 6 



er(xio- 3 ) 

^(xlO- 3 ) 
C^(xl0~ 3 ) 

<P L A°) 



The dependence of the observables B(B — > tttt) and Aqp(B — > tttt) on the new weak phase 
4>%, with definite values of , , and (^ R marked in the legends, is shown in Fig. [TJ We 
can see that both BiB — > tttt) and Acp(B — > tttt) are most sensitive to the Z' contributions 
involving , and Q' , but n °t ^° the ones involving (^ L . Thus, as shown in Figs. UJc) 
and (c'), the Z' effects involving Q' could always be neglected compared with the other ones. 
To be more specific, Our analyses are divided into the following three different cases: 

• Case I: With only C^ R arbitrary, and neglecting all the others; 

• Case II: With and ^ L ' LR arbitrary, but neglecting £j' L ; 

• Case III: Combining with the fitting results from B — > ttK, ttK* and pK decays. 

Corresponding to such three cases, our fitting results for the Z' parameters are summarized in 
Table r2J With these values of Z' parameters as input, our predictions for the observables in 
B — > tttt decays are listed in Table |3l 

Case I: With only R arbitrary, and neglecting all the others. 

In this specific case, as can be seen from Eqs. f|T2|) and ([TBI , the Z' contribution is mainly 
embodied in the QCD penguin sector AC5. From Figs. DJb) and (d), with ~ 0(1O~ 2 ) 
and/or ~ C(10" 2 ) respectively, we can see that B(B° — > tt + tt~) is reduced to be consistent 
with the experimental data at ~ —50°, which is helpful to reconcile the "tttt puzzle". 
However, as shown in Fig. UJb'), the Z' contributions involving a positive ^\ R at 0^ ~ —50° 
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Figure 1: (Color online) The dependence of B(B — > tttt) and Aqp(B — > tttt) on the new weak 
phase 4>\ with the values of £^ L , , Cj" L an d Cd marked by the legends. 
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Table 3: The CP-averaged branching ratios (in units of 10 -6 ) and the CP asymmetries (in unit of 
10 -2 ) of B — > 7T7T decays within the non- universal Z' model with the regulation scheme II for the 
end-point divergence. The theoretical errors of our predictions in Case I, II and III correspond (in 
this order) to the uncertainties of "parameters listed in appendix C and m 5 " and U Z' parameters" . 

Observable Exp. data SM Z' model 

Case I Case II Case III 

B(B~ -> tt-tt ) 5.59±^J 4.52 ±0.70 4.64 ± 0.61 ± 0.0 5.46 ± 0.72 ± 0.83 4.10 ± 0.51 ± 0.12 

B(B° -> 7T+7T-) 5.16 ± 0.22 8.14 ± 1.21 4.52 ± 0.68 ± 0.23 5.71 ± 0.81 ± 0.93 4.66 ± 0.64 ± 0.29 

B(B° -> 7r°7r°) 1.55 ± 0.19 0.95 ± 0.18 1.41 ± 0.20 ± 0.18 1.60 ± 0.23 ± 0.32 1.16 ± 0.16 ± 0.14 

A CP (B~ ->■ tt-tt ) 6 ±5 -0.1 ±0.1 -0.1 ±0.1 ±0.0 0.1 ±0.1 ±0.8 0.0 ± 0.0 ± 0.3 

A C p(B° -»■ 7T+7T-) 38 ±6 26 ±3 32±5±5 29±4±6 30±5±5 

A CP (B° -»■ 7T 7r°) 43lg -45 ±5 -8±4±2 11 ± 3 ± 7 9±3±4 

7T+7T-) -65 ±7 -64 ±7 -53 ±8 ±4 -52±7±4 -50 ± 8 ± 4 

A%$(B° ->• 7T°7r ) 33 ±6 82±8±7 83 ±8± 13 90±5±4 



induce a large negative v4cp(-B° — > 7r 7r ), which is in obvious conflict with the experimental 
measurement 0.43lo24- As a consequence, the case with Z' contributions involving a positive 
at 0^ ~ —50° could be excluded, and we are left with a possible solution to the u im puzzle" 
with only arbitrary as defined by Case I. 

Corresponding to regulation scheme I and scheme II for the end-point divergency respec- 
tively, the final allowed regions for the Z' coupling parameters (^ R and <j)% under the constraints 
from B(B — > tttt) and Aqp{B — > tttt) are shown in Fig. |2j The corresponding numerical results 
are listed in the second column of Table [2J Obviously, one may find that the result of (% R (|^|) 
in scheme II is slightly larger (smaller) than the one in scheme I. 

From the fourth column of Table [3l we can see that most of the observables in B — > -kit decays 
are consistent with the experimental measurements within 2a error-bars. Most importantly, 
our theoretical predictions B7^___ = 1.92 and Rqq = 0.63 agree well with the experimental 
results R%L = 2.02 ± 0.17 and R$ = 0.60 ± 0.08 0, respectively, implying that the Z' 
contributions with (^ R ~ 0.05 and (fij ~ —50° are helpful to resolve the observed "7T7t puzzle". 
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Figure 2: (Color online) The final allowed regions for the Z' coupling parameters Cr and <^ 
in Case I, constrained by the current experimental data on B(B — > tttt) and Acp(B — > tttt). 



It is also noted that, with such parameter values as inputs, our prediction for the mixing- 
induced CP asymmetry A™p(B° — > 7r~7r + )~ —0.53, is in good agreement with the experimental 
measurement —0.65 ± 0.07 [T]. 

In this simplified case, the result Acp{B° — > 7r°7r°) ~ —8% is, however, still ~ 2a lower than 
the experimental data (43^24)%. Even though the Z' contributions with a large R are helpful 
to enhance A C p(B° — > 7r°7r°), such a large value of R is already excluded by B(B° — > tt + ty^) 
as shown in Fig. QJd). 

Case II: With (^ R and ^ L ' LR arbitrary, but neglecting (£ L . 

As discussed in Case I, although most of the observables in B — > tttt decays could be accom- 
modated, it is difficult to reproduce the observed direct CP asymmetry A C p(B° — > 7r°7r°). 
However, as shown in Figs. Q^a') and (b') [], this discrepancy could be possibly compensated by 
the extra Z' contributions involving positive and/or negative with ~ —50°. This 
observation motivates us to consider the second case defined by Case II, where the parameters 
(% R and i^ L ' LR are involved at the same time. 

In this case, the final allowed regions for the four parameters , Cd R and 0^ con- 

strained by BiB — > tttt) and Aqp(B — > tttt) are shown in Fig. [3j and the corresponding nu- 
merical results are listed in the third column of Table [2J Similar to the situation in Case I, we 
1 The plots in the range 0^0^^ 180° could be treated as the ones in —180° ^ (f>^ ^ by shifting the 

,1 • r e LL.LR 1 /-LL,LR 

overall signs of £ d and Q, ' 
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Figure 3: (Color online) The final allowed regions for the Z' coupling parameters £% L , £% R , Cd R 
and <f>j in Case II. The other captions are the same as in Fig. [2J 



again find that the fitting value \<f>^\ in scheme I is a bit larger than the one in scheme II. From 
Fig. El we can see that the value of Cj^ is constrained to be definitely nonzero, which confirms 
our conclusion made in Case I, i.e., the Z' contributions with a NP weak phase 0^ ~ —50° 
and Cj R ~ 0.05 are crucial to reconcile the observed "7T7t puzzle". The constraints on and 
are, on the other hand, not very stringent and both of them could be even equal to zero, 
indicating that they play only a minor role in resolving the u irir puzzle". However, our result 
Acp(B° — )■ 7T 7r°) = 0.11 ± 0.03 ± 0.07 shows that the effects involving ^ L ' LR are still very 
crucial to bridge the large discrepancy of Acp{B° — > 7r°7r°) encountered in Case I. 

As discussed in detail in our previous publication [61], the values of ^ L and £g are found 
to have different signs, which implies that the ratio Ri defined in Eq. ffT"i|) should be negative. 
Since the ratio i?i is independent of the quark flavor d or s, the two parameters ^ L and 
are therefore also expected to have different signs. However, as shown in Figs. 0Ja) anddIV), 
as well as in Fig. [3j a positive (£% R ) is suppressed (even more strongly) by the observalbes 
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B{B~ — > 7r~7r°) and Acp{B° — > ti + tx~)). Thus, after taking into account the fitting results 
of R\ t 3 from Ref. [M], the obtained ranges of Z' parameters (,% L and ^ L,LR could be strongly 
reduced or even be excluded. 

Case III: Combining with the fitting results from B — > nK, %K* and pK decays. 

In our previous work [64J, we have systemically investigated the effects of the family non- 
universal Z' model on penguin-dominated B — > irK, ttK* and pK decays, pursuing possible 
solutions to the observed u ttK puzzle". We found that the model parameter spaces involving 
positive £g L and negative £g (none of them could be neglected) with 0^ ~ —86° are crucial to 
reconcile the puzzle. However, the Z' contributions involving Csi which has a large uncertainty, 
are found to be almost irrelevant. It is interesting to investigate if further information on the 
model parameters could be obtained from the current experimental data on tree-dominated 
B — > Tin decays, which is the motivation of Case III. 

As discussed in Ref. [61], the ratios Ri^ defined by Eq. (fl4|) have already been severely 
constrained, with 

-0.77^ i?i ^ -0.07, -3.7 ^ R 3 < 6.5 . (15) 

Taking these ranges as input, we are then left with only two free parameters 0^ and 

With the current experimental data on B[B — > nn) and AcpiB — > tttt), as well as the ranges 

for i? 13 given by Eq. (1151) as constraints, the final allowed regions for the parameters , , 

and 0^ are shown in Fig. |H The corresponding numerical results for these parameters 

are listed in the last column of Table El As the constraints Eq. (Tl5l) are obtained with the 

regulation scheme II [M], for consistency, in this case we just use regulation scheme II for the 

end-point divergency in our numerical evaluations. 

Comparing Fig. |4] with Fig. [3j one can see that a large part of the parameter spaces of 

^ L,LR allowed in Case II is excluded by the ratios -Ri,3. However, the effects of Ri^ on 

are found to be very tiny. On the other hand, the data on B — > tttt decays could put a 

further constraint on the ratio R3, 4.5 ^ R% ^ 6.5, which implies that the part of parameter 

spaces with ^ R < 4.5 x obtained in Ref. [61] will be excluded. The range of the ratio R± 

2 Of course one could choose any one of £% L , £,^ R and C^ R as the free parameter, and the other two can then 
be reduced from Eq. (Tl4|. It is found that these different choices are irrelevant to our final fitting results. 
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Figure 4: (Color online) The final allowed regions for the Z' coupling parameters ', £% R , 
and cj)^ in Case III. 



remains almost unchanged. As a consequence, the NP parameter spaces are severely reduced 
but still not excluded entirely, which means that both the u wK" and "tttt" puzzles could be 
accommodated simultaneously within such a specific model. 

In order to see the relative strength of the flavour-changing couplings B\ h and B^ b , we define 
the ratio 



s-LL,LR Tr* tjL 

Sd _ \ 

fLL,LR ' I T/* I I Ri 

q,s v td n sb 

-LL,LR 



p — | " ts | | 'db | f-\a\ 

R ds = rr LR ~ 777 75J7 • (IbJ 



With the central values of £ d ' listed in Table El £LL,lr f rom j^ e f jg^j and the CKM param- 
eters listed in Appendix C2, we get numerically 

cLL cLR t /■* 

— = 46 -^— = 4 5 |-*i|=48 (17) 

which implies an interesting relation \B^ b \ ~ I-Bjjj. Thus, in such a family non- universal Z' 
model, the difference of flavour- changing Z' couplings between the quark-level transitions b — > 
dqq and b — > sqq arises only from the new weak phase 0^, with 0^ ~ —50° and 0f ~ —85° [64] . 
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It should be noted that the B® — B d mixing puts a very strong constraint on the b — d — Z' 
coupling B L db [591 EQl EH EH [731 QUI]- As the Z' contributions to B — > -kit decays involve not 
only the b — d — Z' coupling but also the flavor-conserving Z' couplings d — d — Z' (B dd R ) 
and u — u — Z' (B^ R ), one could not extract the constraint on B\ b solely from these decays. 
However, we find that our result (f) d ~ —50° agrees well with the recent combined constraint 
(f)L = -33° ± 45° HE] from B° d - B° d mixing and B s -> n~K + decay. 

Moreover, with the constraint from B Q S —B° S mixing on the b — s — Z' coupling B^ b included, it 
is found that \B db / B^ b \ ~ C(10 _1 ) [S3 EH], which is inconsistent with the relation \B db \ ~ \B^ b \ 
found in this paper. It is noted that the result \B db /B^ b \ ~ (9(10 _1 ) is based on the former 
CDF and DO combined result [102] for the B° s — B° s mixing, however, both the CDF and DO 
collaborations have updated their measurements of the weak phase /3 S very recently, 

{ [0.02,0.52] U [1.08,1.55] CDF [103] , 
0.38l{£jj ± 0.01 DO HD3I, 

which are different from the former CDF and DO combined result [0.27, 0.59] U [0.97, 1.30] |102j . 
Notably the CDF updated measurement [103j agrees with the SM expectation j3 s ~ 0.018 at 
~ la level, while the DO updated measurement [104] agrees with their former combined result. 
As discussed in Ref. [101] . if the lower bound of the CDF updated measurement (3 S ~ 0.02 is 
allowed, one may easily find that the relation |-B^ 6 | ~ \B^ b \ still survives the constraints from 
Bd, s — B<i,s mixings. In such a situation, refined experimental measurements and theoretical 
predictions are therefore expected urgently. 



4 Conclusions 

Motivated by the large discrepancy of the ratio Rqq between theoretical predictions and ex- 
perimental measurements, we have investigated the effect of a family non-universal Z' model 
on the tree-dominated B — > nir decays, pursuing possible resolutions to the observed u irir puz- 
zle" . Moreover, we have also taken into account the fitting results from the penguin-dominated 
B — > irK, irK* and pK decays in Case III, which gives a much stronger constraint on the 
flavour- changing Z' couplings. Our main conclusions are summarized as: 

• The Z' contributions with C^ R ~ 0.05 and <p d ~ —50°, being mainly relevant to the 
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coefficient of QCD penguin operator O5, are crucial to bridge the large discrepancy of 
B(B° — > 7i°7T°) between the theoretical prediction and the experimental measurement. 
The contributions involving the other parameters ^ L,LR and L are almost irrelevant to 
the observed "7T7t puzzle". 

• Combining with the fitting results of Ri^ from the penguin-dominated B — > nK, irK* 
and pK decays, the NP parameter spaces are severely reduced but still not excluded 
entirely. This means that both the "ttK" and u tttt" puzzles could be accommodated 
simultaneously within such a family non-universal Z' model. 

• For all of the three different cases, a new weak phase 0^ associated with the flavour- 
changing Z' coupling B% b , with a value around —50°, is always required for resolving the 
observed discrepancies. 

• The flavour-changing Z' couplings \B% b \ and \B^ b \, corresponding to the quark- level tran- 
sitions b — y dqq and b — > sqq transitions, respectively, are found to be almost equal to 
each other. 

With the upcoming LHC-b and proposed super-B experiments, the data on B-meson decays 
is expected to be more precise |105U106lll07| . which will then severely shrink or totally excluded 
the model parameter spaces. It is also reminded that more refined measurements of mix- 
induced CP asymmetries in these decays are urgently needed to further constrain the Z' coupling 
parameters. 
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Appendix A: Decay amplitudes in the SM with QCDF 



The decay amplitudes for B — > nn decays are recapitulated from Ref. |23j : 

3 



V2A™ 



A 



SM 
B°- 



Y V Pb V pd A 7Tn $pu {Oil + «2) + 2 («3,EW + a 

)=Ji,C 

V PbV p *d A nn [S p u («i + Pi) + a% + a* EW + (3^ + 2(3, 

(/^3,EW — ^4,Ew) 



(19) 



p=u,c 

1 

~2 



(20) 



aSM 



r 3 1 

Yl V Pb V P *d A ™ $pu («2 - 0l) - «4 + a«3,EW + a«4,EW 



p=u,c 

1 



+ o (^3,EW /^4,EW; 



(21) 



where explicit expressions for the effective coefficients af = a p i {MiM 2 ) and (3f = /3f(M!M 2 ) 
could also be found in Ref. [22] • It should be noted that the hard-spectator terms Hi ap- 
pearing in of and the weak annihilation terms A* appearing in /3J should be replaced by our 
recalculated ones listed in Appendix B. 



Appendix B: The hard-spectator and annihilation kernels 
with the infrared finite gluon propagator 

With the infrared finite gluon propagator to cure the end-point divergence, the hard-spectator 
kernels in a general B — > PP decay can be expressed as [76] 

w\ b mim 2 f 1 , , Ac ^s{g 2 ) ^ , v r $m 1 (y) Ml <t> mi (y) ] 

H ' (MlM2) = A^-Jo dxd ^— *"«)*"-( I Hx( S + ^)/0 +r x x(y + ^)/i)\' 

(22) 



for the insertion of operators Qi=i-4, 



9,10, 



Am 1 m 2 Jo £ 



Mi 



<t>mi{y) 



x(y + uJ 2 (q 2 )/0 x x(y + u 2 (q 2 )/O i 

(23) 
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for Qi=5.7, and Hi(MiM 2 ) = for Qi=e,s- When both Mi and M 2 are pseudo-scalars, the final 
building blocks for annihilation contributions can be expressed as [76] 



A\ = 7r / dxdya s (q 2 
Jo 



X 



+ 



xy — u 2 (q 2 ) + ie 



A2 = n dxdya s (q 
Jo 



(xy — u 2 (q 2 ) + ie) (1 — xy) (xy — u 2 (q 2 ) + ie)x 
rfrf cf> mi (y)<f> m2 (x)} > 

1 



— w 2 (g 2 ) + ie 



(xy — ui 2 (q 2 ) + ie)(l — xy) (xy — uj 2 (q 2 ) + ie)y 
2y 



$ Ml (y)^M 2 (a;) 

(24) 

<5> Ml (y)®M 2 (x) 
(25) 



4 = Try d X dya s ^ UTy _ u2{q2)+ . e)ii _ x _ y; 



-rf mi (y)$ A/2 (x) 



2x 



rf 2 (x)</> m2 (x)$ Ml (y)|> , 



(xy — u 2 (q 2 ) + ie)(l — xy) x 



(26) 



A{ = 4 = 0, 



2(1 + x) 



7T / dxdya s (q 2 ) ] 
Jo [(xy — UJ z (q z ) + te)x ' 



r^cf> mi (y)^ M2 (x) 



2{ l^ V \ . ; rf (x)0 m2 (x)$ Ml (y) 



(xy — uj 2 (q 2 ) + ze)y x 



(27) 



(28) 



Appendix C: Theoretical input parameters 



CI. The numerical results of Wilson coefficients 

The numerical results of Wilson coefficients in the naive dimensional regularization (NDR) 



scheme at the scale fi = m b ([i h = ^Jk h m h ) are listed in Table SJ For simplicity, we have defined 



x = \VtbV t *d\ C LL i<f> L y = \VtbV*d\ C LR i<f> L 



y/ _ \VtbVtd\ fLL %<j>L 

v tb v td Q 



V l _ \VtbVtd\ r LR i<p L 

v tb v td Q 



(29) 



The values at the scale with m b = 4.79 GeV and = 500 MeV, should be used in the 
calculation of hard-spectator and annihilation contributions. 



19 



Table 4: The Wilson coefficients Cj within the SM and the Z' model in NDR scheme at the scale 
/i = nib and p,h = sj A^mf,, respectively. 



Wilson 


fj, = m b 


Ph = V^hm b 


coefficients 


C 


i 




AC? 




AC? 


Ci 


1 


075 




-0.006X 


1.166 


-0.008X 


c 2 


-0 


170 




-0.009X 


-0.336 


-0.014X 


c 3 





013 


0.05X - 0.01Y - 2.20X' - 


- o.osy' 


0.025 


0.11X - 0.02Y - 2.37X' - 0.12Y' 


c 4 


-0 


033 


-0.13X + 0.01Y + 0.55X' - 


h 0.02Y' 


-0.057 


-0.24X + 0.02Y + 0.92X' + 0.09F' 


c 5 





008 


0.03X + 0.01Y - 0.06X' - 


- 1.83Y' 


0.011 


0.03X + 0.02Y - 0.10X' + 0.09F' 


c 6 


-0 


038 


-0.15X + 0.01Y + 0.1X' 


- 0.6Y' 


-0.076 


-0.32X + 0.04Y + 0.16X' - 1.26Y' 




-0 


015 


4.18X-473Y + 0.25X'- 


h 1.27Y' 


-0.034 


5.7X - 459Y + 0AX' + 1.7Y' 


C%ja em 





045 


1.18X - 166Y + 0.01X' + 0.56Y' 


0.089 


3.2X - 355Y + 0.2X' + 1.5Y' 


Cg/a em 


-1 


119 


-561X + 4.52Y - 0.8X' 


+ 0.4Y' 


-1.228 


-611X + 6.7Y - 1.2X' + 0.6Y 7 


Cio/Oem 





190 


118X - 0.5Y + 0.2X' - 


- 0.05F' 


0.356 


207X - IAY + 0.5-X 7 - 0.1Y' 




-0 


297 






0.360 




C 8g 


-0 


143 






-0.168 





C2. CKM matrix elements 

For the CKM matrix elements, we adopt the Wolfenstein parameterization |1U8] and choose 
the four parameters A, \, p and 77 as [109] 

A = 0.7981°;°??, A = 0.225211™, p = 01411°;°!?, V = 0.340 ± 0.016, (30) 
with p = p (1 — and fj — rj (1 — ^-). 

C3. Quark masses and lifetimes 

As for the quark masses, there are two different classes appearing in our calculation. One type 
is the current quark mass which appears in the factor r^ 1 through the equation of motion for 
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quarks. This type of quark masses is scale dependent and denoted by m q . Here we take 

m s (fi)/m q (fi) = 27.4±0.4 [TJD], m s (2GeV) = 87±6MeV [110], m b (m b ) = 4.20+°;^ GeV 

(31) 

where WT g (/i) = (fn u + m d )(/i)/2, and the difference between u and d quark is not distinguished. 

The other one is the pole quark mass appearing in the evaluation of penguin loop diagrams, 
and denoted by m q . In this paper, we take 



m u = md = m s = 0, m c = 1.61+^2 GeV, m b = 4.79+^ GeV. 



1+0.19, 



(32) 



As for the B-meson lifetimes, we take [12] tb u = 1.638 ps and tb c1 = 1.530 ps, respectively. 



C4. Input parameters related to mesons. 

In this paper, we take the heavy-to-light transition form factors and the decay constants as 



F,f^(0) = 0.258±0.031 [mj E2J E3], f B = (216±22) MeV [EH! , f n = (130.4±0.2) MeV 

(33) 



The light-cone projector operator of light pseudo-scalars in momentum space reads [23 ] 1115] 

(34) 

where fp is the decay constants, and \ip = m b r^/2 with the chirally-enhanced factor defined 



^75 $p(x) - /ip7 5 2 I (j) p {x) 



as 



w _ 2ml 

x m 6 (/i) 2m ? (/i) ' 



(35) 



where the quark masses are all running masses defined in the MS scheme. For the light-cone 
distribution amplitude of light pseudo-scalars, we use their asymptotic forms |116[ I117[ I118[ 1119] 



§p(x) = 6x(l - x) 



'x) = l. 



(36) 



As for the B-meson wave function, we take the form |120j 



MO = N B £(1 - £)exp 



M, 



B 



(37) 



-M B - m b 

where £b = 1 — m b /MB, and Nb is the normalization constant to insure that J* d^s(0 = 1- 
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